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Abstract. This paper presents recent results obtained by the authors (partly in collabo- 
ration with A. Dabrowska) concerning expansions of zonal functions on Euclidean spheres 
into spherical harmonics and some applications of such expansions for problems involving 
Fourier transforms of functions with rotational symmetry. The method used to derive the 
expansion formula is based entirely on differential methods and completely avoids the use 
of various integral identities commonly used in this context. Some new identities for the 
Fourier transform are derived and as a byproduct seemingly new recurrence relations for the 
, classical Bcsscl functions arc obtained. 

43 ■ K e U words: spherical harmonics; zonal harmonic polynomials; Fourier-Laplace expansions; 

special orthogonal group; Bcssel functions; Fourier transform; Bochncr identity 



2000 Mathematics Subject Classification: 33C55; 42B10; 33C80; 44A15; 44A20 



This paper is dedicated to Professor Jacques Faraut 
on the occasion of his 65-th anniversary. 



43 
+-> 

a 

S 

> 

•i-H 

X 

1 Introduction and preliminaries on spherical harmonics 

Many parts of analysis on Euclidean spaces and in particular the theory of spherical harmonics 
provide an elegant and instructive application of group theoretical concepts to various questions 
of classical function theory. It suffices to mention points like: coincidence of eigenspaces of 
spherical Laplacian with irreducible components with respect to the natural action of the group 
SO(d) within the I? space on the sphere; the role of classical orthogonal polynomials, i.e. 
Gegenbauer polynomials, as reproducing kernels for the spaces of spherical harmonics of a given 
degree, or more generally, as providing an explicit construction of symmetry adapted basis 
functions for those spaces (cf . ^2] ) > ° r the connection of the Fourier transform on the Euclidean 
space to the Hankel transform obtained via restriction to SO(cZ)-finite functions and various 
integral identities of the Hecke-Bochner type resulting there from. 

The goal of the present paper is to present under one cover recent developments within 
this circle of ideas, obtained in |1 11 HI El by the present authors, partly in collaboration with 
A. Dabrowska. The central point of this paper is a novel form of spherical expansion of zonal 
functions on Euclidean spheres which we derive by purely differential methods. We show how it 
implies certain general formulae for the Fourier transform of SO(d)-finite functions, including 



2 



A. Bezubik and A. Strasburger 



recent generalizations of the famous Bochner formula and also derive certain function theoretic 
consequences. 

1.1 Preliminaries and notations 

We denote by (x \ y) the standard (Euclidean) inner product of points x, y E W 1 and by r = 
|x| = (x | x) 1 / 2 the corresponding length (or radius) function. The unit sphere in M. d is denoted 
by S d ~ l and for x 7^ we shall often write x = r£ with £ E S* ' -1 . We shall assume d > 3 and 
frequently use a constant a defined by a = (d — 2) /2. The sphere is regarded as a homogeneous 
space of the group K = SO(d) of (proper) rotations in M rf and for a point i] E S^" 1 its isotropy 

group -K^ C K can be identified with the group SO(d — 1). 

d 

Let A = ^ d 2 /dx 2 denote the Laplacian in R d . We let V 1 = V l (M. d ) denote the space 
i=i 

consisting of homogeneous polynomials of degree I on R rf and denote the kernel of A in V 1 
by H l . This latter space consists of harmonic and homogeneous polynomials of degree I and its 
dimension is given by 



dim W 



, _ 2{l + a)T{2a + l) 
~ r(Z + l)r(2a + l) 



where T(z) denotes the Euler gamma function. Both these spaces are invariant under the natural 
action of the group K = SO(d) on functions on M. d given by k ■ f(x) = f(k~ 1 x) for k £ K, 
x £ M. d . It is known that this group acts irreducibly in TL l for each /, and the representations so 
obtained are inequivalent for I 7^ I'. 

The surface spherical harmonics of order I are by definition the restrictions of elements 
from TL l to the unit sphere S^ -1 , and since the restriction map commutes with the action of 
rotations the spaces of surface spherical harmonics of any fixed order are invariant and irreducible 
under K. 

We endow V 1 = V l (M. d ) with an inner product defined by the formula 



(p\Q) := / P(OQ(OMO, (i) 

where da(-) is the Euclidean surface measure on the unit sphere 5 d_1 normalized so that the 
total area of the sphere is 1. In fact, the integral on the right hand side extends to the natural 
inner product on L 2 (S d ~ l , da) and it is known that the spaces of spherical surface harmonics of 
different orders are orthogonal to each other with respect to this inner product. 

We recall that a function /, defined on the unit sphere S d ~ l , is said to be a zonal function 
(relative to a point n £ S d ~ l ) if it is invariant with respect to the isotropy group K v of rj. Any 
such function depends in fact on the inner product (£ | rf) only, and so there exists a function <j) 
defined on the closed unit interval [—1, 1] C M so that 

f(0=m\v)), for all£ G S d ~\ (2) 
We shall call <f> the profile function of /. 



1.2 The canonical decomposition of homogeneous polynomials 

It is well known that every homogeneous polynomial of degree / can be represented as a sum of 
products of harmonic homogeneous polynomials with even powers of the radius, cf. equation © 
below, and this decomposition is usually called the canonical decomposition of homogeneous 
polynomials. While the general form of the canonical decomposition, i.e. equation (j2J), is stated 
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in all sources concerned with this matter, the explicit formula for harmonic components of 
a given homogeneous polynomial is seldom quoted. The only source known to us, where the 
formula can be found, is the paper [7j of Lucquiaud. We record it here since it is essential for 
considerations to follow, in particular, in establishing formula © giving expansion of elementary 
zonal polynomials of the form (x \ rj)\ where rj G S 1 is an arbitrary unit vector. 

Theorem 1 (The canonical decomposition). The space V 1 decomposes orthogonally as the 
sum V 1 = @\^r 2k TL l ~ 2k and the decomposition is invariant with respect to the action of the 
group SO(d). Explicitly, every polynomial P G V 1 may be written as 

[1/2] 

P = J2r 2k h^ 2k (P), where fc,_ 2fc (P) G H l ~ 2k (3) 

k=0 

are called the harmonic components of P and are given by 

[l/2]-k 

h-2k{P)= Yl ^{k)r 2j ^ j P (4) 

3=0 

with the coefficients e l j(k) determined by 

, j (a + l-2k)T(a + l-2k-j) 

e i {k) ~ [ ij &+jk\j\T(a + l + l-k) ' 1 j 

The maps P i— > r 2k hi_2k(P) G P l are projections onto SO '(d) -irreducible subspaces of V 1 com- 
muting with the action of the group SO(d). 

2 Expansions of zonal functions 

2.1 Elementary zonal functions and reproducing kernels 

To proceed, we need to recall an explicit formula for the Gegenbauer polynomial Cf of degree I 
and index a, given for example in jSJ Section 5.3.2] or |131 Chapter 9.2], which reads as follows: 

Applying the formulae for the canonical decomposition given in © and © to elementary 
zonal polynomials (x \ rj) 1 , where r\ G S^" 1 and / is a nonnegative integer, we obtain 

(x | n) 1 = 2- l T(a)T(l + l)\x\ l Y, kT ( a + l + 1 _ k) C ?-2k((Z 1 17)) , x= \x\£. (6) 

In particular, the spherical harmonic obtained by restricting to the unit sphere the harmonic 
component hi(P v ) of the highest degree of -Pjj(x) = (x \ rf) 1 , which is given by the formula 

Ji,(P„)(fl = 2- l r(a)T(l + l )f ^±iL_CT((e I V)), 

plays an important role in the group representation theoretic interpretation of the decomposi- 
tion (|TJ). With normalization given by 



z l v (0 = [c?(iT 1 cm\v)) 
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it satisfies 

dimW' / Z l (OP(Oda(0=P(v)- 

Because of this property, is called the reproducing kernel for the space 7i l . Moreover, it 

is uniquely (up to a scalar multiple) determined by the property of being invariant under the 
action of the isotropy subgroup K v of the point rj 6 S d_1 . 



2.2 A differential formula for expansions of smooth zonal functions 

Now recall that every square integrable function on the sphere can be written as a series 
of spherical harmonics (the Fourier-Laplace expansion). For zonal functions this expansion, 
thanks to an easy group representation theoretic argument, reduces to 

f(0 = E fraZ'viO, where f m = dim H m / /(£)^(0 MO- 

m=0 ^ Sd 1 

By taking into account the invariance under of the integrand and equation (J2J), the coefficients 
may be expressed as integrals 

what reduces the problem of spherical expansion of zonal functions to the expansion of profile 
functions <f> with respect to the (orthogonal) system of Gegenbauer polynomials. The following 
result shows that the coefficients f m of the expansion can also be expressed in terms of the 
coefficients of the Taylor expansion of the profile function <f>, provided the latter satisfies suitable 
regularity assumptions. 

Theorem 2. Assume <f) : [—1,1] — > C has the Taylor expansion 4>{t) = £ ^-^t m that is 

absolutely convergent on the closed interval [—1, 1], and let /(£) = </>((£ | r])) be the zonal function 
on the sphere S^" 1 corresponding to 4>. Then the spherical Fourier expansion of f{£) is given by 

oo 

/(0 = T(a + / m dimW m Z-(0, (7) 

m=0 

where the coefficients of the expansion can be expressed as 

Jm- 2^ 2 m + 2k k\T(a + m + k + 1) ' 

A detailed proof of this result is contained in the forthcoming paper [2] of the authors, and 
here we shall present its main line only. It consists in substituting the expansion formula © into 
the Taylor series of <f> and rearanging terms so that to group together the terms corresponding 
to the Gegenbauer polynomials of a given degree. This procedure requires some estimates on 
the coefficients which assure the absolute convergence of the double series, which are presented 
in detail in [3]. 

Below we briefly present two applications of this expansion for obtaining new derivations of 
some classical results. 
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2.3 The plane wave expansion 

A well known and very useful instance of the expansion (J7J) is the so called plane wave expan- 
sion, giving a representation of the exponential function e l ^ x ^ as a series of zonal harmonic 
polynomials. The expansion involves the Bessel functions of the first kind of order v S C with 
Rez^ > — 1, which are given as 



For some purposes the results are better expressed with the aid of the so called spherical Bessel 
functions, which are given as 

j„(t) = r(v + i)(^\ " j v (t). 

We point out that all classical proofs of this formula known to us are obtained by applying 
certain integral identities of the Hecke-Bochner type, as in With the use of (J7| we get it 

by direct differentiation of the exponential e irt and comparison with the power series expansion 
of the Bessel function (|8|). 

Corollary 1. For an arbitrary unit vector n € S^ 1 C M. d and x £ M. d , the plane wave e^ x '^ 
admits the following expansion 

e i(*k) = dim ft" l " + j ( r -) m j a+m (r)Z?(Q, x = ri, (G^ 1 . 

^-^ 1 (a + m + 1 \2/ ' 

m=0 v ' 

T/ie series converges absolutely on each sphere of radius r and uniformly with respect to £, 77 G 

s d -\ 



2.4 The generating function of Gegenbauer polynomials 

Take a fixed element i ^ from the unit ball in W 1 and write x = rr), rj 6 S d ~ l . The function 

S^" 1 3 £ h- |f - x|" 2Q = (1 - 2r(£ | r?) + r 2 )"° 

is clearly a zonal function with pole at i] corresponding to the profile function (1 — 2rt + r 2 )~ Q . 
The expansion resulting from (JJJ) has the form 

1 {U) } ^ r(2a)r(m + i) 

and can be reduced to the familiar formula for the generating function of Gegenbauer polynomials 

OO 

(l-2rt + r 2 y a =Y,r m C«(t), 

m=0 



as given e.g. in IB]- 
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3 The Fourier transforms of SO(d)-finite functions 

In this section we shall present some applications of the plane wave expansion given in Corollary^ 
to problems of Fourier analysis in Euclidean space. The subject is related to the well known 
Bochner formula (c.f. e.g. 0113), which describes the restriction of the Fourier transform to the 
space of functions of the form f(r)P(x), where / is a radial function and P G T~i l a homogeneous 
harmonic polynomial. The result is again a product of the same harmonic polynomial P with 
a radial function, which is expressed as the so called Hankel transform of the original radial 
factor and is given by the formula 1)12(1 below. The result we present in Corollary 0] generalizes 
that relation for the case of arbitrary homogeneous polynomials and is a combination of results 
obtained by the authors in [111 12] and by F.J. Gonzalez- Vieli in [H]. 

3.1 The case of SO(d)-finite functions on the sphere 

We take the Fourier transform of suitable regular (e.g. L 1 (IR c( )) functions on M. d as defined by 



Observe that the above definition of the Fourier transform also makes sense in the case when 
/ is a function defined on the unit sphere — in this case it can be regarded as the Fourier 
transform of a measure supported on the sphere. Especially interesting is the case, when the 
measure on the sphere comes from restricting a homogeneous polynomial to the sphere, since 
this way we obtain a so-called SO(d)-finite measure (since its SO(<i)-translates span a finite 
dimensional subspace) . 

Theorem 3. If P G V , then the Fourier transform of the measure P(£)da(£) with support on 
the unit sphere S"^" 1 is given by the following equivalent formulae 



with h[-2k(P) denoting the harmonic components of P as in equation ©. 

The formula (|1U|) has been derived by F.J. Gonzalez |S] and the equivalence of those two 
forms was observed by the present authors and A. Dabrowska in [2]. 

Theorem |3] implies the following interesting function theoretical corollary. By comparing 
the two expressions © and (|10[) from this theorem one may derive the following multi-step 
recurrence relation for spherical Bessel functions. 

Corollary 2. If a > is a half odd integer (or an integer) then for any I G Z + the following 
relations hold 





(9) 




(10) 
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for s = 1 



...,[1/2]. 
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In terms of the Bessel functions of the first kind, this is expressed by 

1 ( 2 Y t i \ - ST T(a + / - fc - sjTja + I + 1 - 2k) , f . 

J.{r) J ^ l ~^ ~ k \(s - k)\T(a + 1 + 1- k)T(a + I - 2k) J ^~^ W 



It is interesting to note that these latter relations unify several classical relations satisfied by 
Bessel functions like ^ equation (4.6.11)] or [SJ Chapter 3.2.2]. In fact, taking into account the 
relation 

1/2 



J-1/2W) = ( 2 71 "* J cost 



which follows directly from the expansion (JHJ), and the familiar recurrence relations satisfied 
by J u (t), namely 

{t- v J v {t)) = {-i)h-»- l j v+l {t) 

one can derive from 1)11(1 the following relations. 

Corollary 3 (Finite expansions of Bessel functions). The Bessel functions of integer order 
satisfy the relation 

1 /2\ n 



f2\ n 1 

[t) J ^ = ^ {n+mn - k) M ^ where ek 



n! V*/ {n + k)\{n-k)\ z " w ' 1 2, fc = l,2,3..., 

u;/izZe £/iose o/ /ia// ode/ integer order satisfy 



■4+1/2 CO = \/— * 




7TX 



[n/2] 

sin(i — nn/2) 



(-l) fc (n + 2A;)! 
Q (2fc)!(n-2fc)!(2t) 2fc 



[ ( ^ /2] (-!)*(„ + 2* + 1)! 
+ cos(i - nvr/2) ^ ( 2 fc + l)!( n - 2fe - l)!(2t)^+ i f 

3.2 The generalized Bochner formula and the Hankel transform 

For our last topic recall that for suitable functions on the nonnegative real axis M+, say, with <p 
belonging to the Schwartz space <S(R), one defines the Hankel transform by the formula 

n — u roo 

H»(m) = tv , T\ / 0(s)j„(st)s 2 » +1 ds. (12) 



r(i/ + i) 7 

Theorem above now immediately implies that the generalized Bochner formula which was 
previously obtained in can also be expressed by a pair of equivalent formulae. 

Corollary 4 (Generalized Bochner identity). If P £ i/ien i/ie Fourier transform of the 
function f(\x\)P(x) is given by the following expressions: 

, [»/2] 

(2vr)-i / /(|x|)P(x)e^)dx = i 1 V(-l) fc fl-a+l-2fc(s 2fc /)(|y|)/»l-2fc(i')(y) 



k=0 

V/2] (-l) k 

•'E Wr^- fc (/)(bl)A fe p( y ) 

fc=0 
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This implies in turn the following: 

Corollary 5 (Periodicity relation for the Hankel transform). For any <j) & S(M.) with a 
and I satisfying conditions of Corollary |2] above, the Hankel transform satisfies the following 
relation 

H a+l {cp){t) =2(a + l- 1)F Q+J _i(0)(i) - H a+l _ 2 (s 2 4>) (t). 
4 Conclusion 

In this paper we have demonstrated that differential identities for homogeneous polynomials 
like those implied by equations can be effectively used for solving problems in harmonic 

analysis, which so far have been approached by means of integral identities of the type of Hecke- 
Bochner formula. In our opinion other possibilities of using that approach should certainly be 
further explored. 
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